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ABSTRACT 

We derive rotation curves for four nearby, fow-mass spiral galaxies and use them to constrain the 
shapes of their dark matter density profiles. This analysis is based on high-resolution two-dimensional 
Ha velocity fields of NGC 4605, NGC 5949, NGC 5963, and NGC 6689 and CO velocity fields of 
NGC 4605 and NGC 5963. In combination with our previous study of NGC 2976, the full sample of 
five galaxies contains density profiles that span the range from Odm = to = 1.20, where aum is 
the power law index describing the central density profile. The scatter in a^M from galaxy to galaxy 
is 0.44, three times as large as in Cold Dark Matter (CDM) simulations, and the mean density profile 
slope is ckoM — 0.73, shallower than that predicted by the simulations. These results call into question 
the hypothesis that all galaxies share a universal dark matter density profile. We show that one of 
the galaxies in our sample, NGC 5963, has a cuspy density profile that closely resembles those seen in 
CDM simulations, demonstrating that while galaxies with the steep central density cusps predicted 
by CDM do exist, they are in the minority. In spite of these differences between observations and 
simulations, the relatively cuspy density profiles we find do not suggest that this problem represents a 
crisis for CDM. Improving the resolution of the simulations and incorporating additional physics may 
resolve the remaining discrepancies. We also find that four of the galaxies contain detectable radial 
motions in the plane of the galaxy. We investigate the hypothesis that these motions are caused by 
a triaxial dark matter halo, and place lower limits on the ellipticity of the orbits in the plane of the 
disk of 0.043-0.175. 

Subject headings: dark matter — galaxies: dwarf — galaxies: halos — galaxies: individual (NGC 2976; 

NGC 4605; NGC 5949; NGC 5963; NGC 6689) — galaxies: kinematics and dynam- 
ics — galaxies: spiral 



1. INTRODUCTION 

Over the last several years, one of the most persis- 
tent problems confronting the Cold Dark Matter (CDM) 
model has been the dichotomy between observed galaxy 
density profiles and those seen in cosmological simula- 
tions. The simulations generate dark matter halos with 
central density cusps of p cx or steeper (e.g., Navarro, 
Frenk, & White 1996; Moore et al. 1999). Observa- 
tions of dwarf and low surface-brightness (LSB) galaxies, 
though, have usually shown that density profiles with 
shallow central cores, where the density is nearly con- 
stant with radius, fit the data better than cuspy profiles 
(Moore 1994; Burkert 1995; Blais-Ouellette et al. 1999; 
Blais-Ouellette, Amram, & Carignan 2001; de Blok et al. 
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2001a; Borriello & Salucci 2001; de Blok, McGaugh, & 
Rubin 2001b; de Blok & Bosma 2002; Salucci, Walter, & 
Borriello 2003; Weldrake, de Blok, & Walter 2003; Simon 
et al. 2003). This disparity is important because den- 
sity profiles represent some of the strongest available con- 
straints on the CDM simulations at small spatial scales. 

Because most observational errors tend to make den- 
sity profiles look shallower than they actually are, 
whether this disagreement is real or only apparent has 
been a point of contention. Recently, van den Bosch 
et al. (2000) and van den Bosch & Swaters (2001) put 
forth the argument that most existing rotation curves 
have neither the spatial resolution nor the velocity reso- 
lution necessary to test the simulations adequately. Swa- 
ters et al. (2003) extended this argument with simu- 
lations showing that even high-resolution data may be 
subject to numerous systematic effects that can make 
density profiles appear artificially shallow. However, de 
Blok & Bosma (2002) and de Blok, Bosma, & McGaugh 
(2003) use observations and simulations to reach the op- 
posite conclusion: that systematic uncertainties do not 
significantly affect most observed rotation curves. In an 
attempt to resolve this debate, we began an effort to 
measure the density profiles of a sample of nearby dark 
matter-dominated galaxies at very high spatial and ve- 
locity resolution, using improved techniques to minimize 
the importance of systematic uncertainties. 
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In two previous papers, we reported on our results for 
tlie dwarf spiral galaxies NGC 4605 (Bolatto et al. 2002) 
and NGC 2976 (Simon et al. 2003). We found that 
NGC 2976 contains a nearly constant density core, while 
NGC 4605 has a density profile that is intermediate be- 
tween the usually observed constant density cores and 
the cusps predicted by simulations. In this paper, we 
present similar analyses of three more nearby galaxies, 
NGC 5949, NGC 5963, and NGC 6689. Since our earher 
study of NGC 4605 did not include a two-dimensional Ha 
velocity field, which prevented us from confirming the ex- 
istence of the radial motions suggested by the CO data, 
we also update that work with an Ha velocity field, addi- 
tional CO mapping, and new optical imaging. As before, 
we use high-resolution two-dimensional velocity fields to 
constrain the overall mass distribution, and optical and 
near-infrared imaging to model and remove the baryonic 
contribution to the rotation curve. 

In the following section, we describe our observations 
and data reduction. We briefiy explain our methods of 
modeling stellar disks and constructing tilted-ring mod- 
els of the velocity fields, and refer the reader to Simon et 
al. (2003) for additional details. In §3, we present the 
results of our rotation curve analysis and fit the rotation 
curves with various functional forms. In §4 we discuss 
the most important implications of this work. We con- 
sider the effect on the derived density profiles of relaxing 
the assumption that dark matter halos are spherical in 
§5 and we present our conclusions in §6. 

2. OBSERVATIONS, DATA REDUCTION, AND 
METHODOLOGY 

The primary objective of this study is to derive very 
accurate rotation curves (and hence density profiles) of 
low-mass,'^ dark matter- dominated spiral galaxies. We 
use two-dimensional velocity fields to measure the rota- 
tion curves in order to avoid the systematic problems as- 
sociated with long-slit spectroscopy (Swaters ct al. 2003; 
Simon et al. 2003). Where possible, we obtain velocity 
fields at multiple wavelengths (e.g., Ha and CO) as a 
further guard against systematic errors. Modeling and 
removing the stellar disk is an important step in the 
derivation of the dark matter density profile, so we use 
multicolor imaging to obtain the best available estimate 
of the stellar mass-to-light ratio. 

Our target galaxies were originally selected from the 
CO survey of nearby dwarf galaxies by Leroy et al. (in 
preparation) in the hope that we could map the CO emis- 
sion in the centers of the galaxies with the BIMA inter- 
ferometer. Only three galaxies in the survey (NGC 2976, 
NGC 4605, and NGC 5963) contained enough CO to pro- 
duce useful interferometer maps, so we also added two 
more galaxies that were similar in mass (with rotation 
velocities of approximately 100 km s^^), distance (ap- 
proximately 10 Mpc), inclination, and luminosity to the 
CO-rich dwarfs. 

2.1. Target Galaxies and Observations 

^ Wc refrain from using the term dwarf galaxy here to avoid 
confusion caused by the range of definitions of the term in the 
hterature. We instead refer to our targets as low-mass galaxies, by 
which we mean that they are signifieantly less massive than giant 
galaxies like the Milky Way and thus can be expected to be dark 
matter-dominated. 



The properties of NGC 4605 arc described by Bolatto 
et al. (2002). This late type galaxy has an absolute mag- 
nitude of Mb = -17.7 and it is located 4.26 ± 0.64 Mpc 
away (M. Pierce 2001, private communication). While 
NGC 4605 is classified as a barred galaxy (SBc), detailed 
images do not reveal evidence for the presence of a bar. 

NGC 5949 is an Sbc galaxy of similar luminosity 
{M-B = —18.2) at a larger distance. Since this galaxy was 
included in the Spiral Field I-band (SFI) Tully-Fisher 
survey, we can use the observed parameters given by 
Haynes et al. (1999a, b), corrected for Galactic and in- 
ternal extinction and for turbulent broadening of the H I 
line, and the Tully-Fisher Relation determined by Gio- 
vanelli et al. (1997) to calculate a distance of 14.0 ± 2.4 
Mpc. 

NGC 5963 is an Sc galaxy associated with the 
NGC 5866 group (Fouque et al. 1992). The nearest large 
galaxy in the group is NGC 5907, at a projected distance 
of 430 kpc, so it is unlikely that NGC 5963 is currently 
interacting with its neighbors. It has a heliocentric re- 
cession velocity of 654 ±3.1 km s~^. There are no pho- 
tometric distance determinations for this galaxy in the 
literature, so it is reasonable to assume that it lies at the 
distance of 13 AIpc implied by the Hubble fiow"* after cor- 
recting for the motion of the Local Group towards Virgo. 
The uncertainty on this distance is probably of order 3 
Mpc. At 13 Mpc, NGC 5963 has an absolute magnitude 
of Mb = -17.8. 

NGC 6689 (also called NGC 6690) is fisted as SBcd 
in NED and SBc in LEDA, but it is classified as an un- 
barred Sc galaxy in the UGC (Nilson 1973), and Sd in 
the Third Reference Catalog of Bright Galaxies (de Vau- 
couleurs et al. 1991, hereafter RC3). No bar is evident in 
any of our optical or near-infrared images. We conclude 
that the confiicting classifications are due to the rela- 
tively high inclination of the galaxy and that NGC 6689 
does not contain a bar. The galaxy is rather similar 
to the others (but more inclined), with a distance of 11 
Mpc (again using the Virgocentric-fiow-corrected veloc- 
ity) and an absolute magnitude of Mb = —17.6. 

2.1.1. Ha Observations and Reductions 

Our Ha observations were obtained on 2002 March 
20-21, 2002 May 25-28, and 2003 April 15 at the 3.5 

m WIYN telescope with the DensePak fiber array. See 
Harden, Sawyer, & Honeycutt (1998) and Simon et al. 
(2003) for details about the instrument and spectrograph 
setup, respectively. 

We observed the galaxies at five to eleven positions, 
depending on the spatial extent of their Ha emission. 
These observations resulted in 350-650 independent ve- 
locity measurements for each galaxy. Exposure times 
ranged from 10 minutes to 60 minutes per position. The 
angular resolution of these data is 4" and the velocity 
resolution is 13 km . 

The DensePak data were reduced as described by Si- 
mon et al. (2003). The only significant change in 
the data reduction resulted from the replacement of the 
CuAr comparison lamp with a ThAr lamp in 2002 May. 
The very bright Th lines near the wavelength of Ha al- 
lowed us to improve the accuracy of the wavelength so- 

^ We use the Hubble Space Telescope Key Project value of Ho = 
72 km Mpc~^ (Preedman et al. 2001) for the Hubble constant. 
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lution by about a factor of two (to ^0.2 km s 

2.1.2. CO Observations and Reductions 

Our (J = 1 0) observations of NGC 4605 

and NGC 5963 were acquired using the B, C, and D 

configurations of the 10-clcnicnt BIMA array (Welch et 
al. 1996) between April 2001 and March 2002. The total 
integration time for each galaxy was ^ 80 hours, much of 
which was in the most extended (B) configuration. The 
CO emission in NGC 4605 extends beyond the BIMA 
primary beam diameter of ~ 100", so we constructed a 
mosaic of observations made at several positions along 
the major axis of the galaxy. (The CO observations of 
the central field used in this paper are the same as those 
of Bolatto et al. 2002, and the other major axis fields are 
from the new mosaic.) Because of the greater distance 
of NGC 5963, the CO emission in that galaxy is much 
more compact and fits easily inside the primary beam. 
Observational setup and data reduction were identical 
to that described by Simon et al. (2003). Beamsizes 
and sensitivities were 5'.'8 x 5'.'1 (120 x 105 pc) and 24 
mJy beam"^ for the NGC 4605 central field, and 15'.'1 x 
13'.'4 (312 X 277 pc) and 70 mJy beam"^ for the outlying 
fields of the mosaic. For NGC 5963, the beamsize was 
5'.'8 X 5'.'3 (370 x 331 pc) and the sensitivity was 31 mJy 
beam-i. We detected CO in NGC 5949 and NGC 6689 
with observations at the UASO 12 m telescope (Leroy 
et al., in preparation), but the emission was not bright 
enough to map with BIMA. 

2.1.3. Optical and, Near-IR Imaging and Reductions 

We observed NGC 4605 and NGC 5963 with B, V, R, 
and I filters at the 1.8 m Perkins Telescope at Lowell 
Observatory on 2002 February 11. NGC 5963 was also 
observed with much longer exposure times in B and R 
with the 1 m Nickel Telescope at Lick Observatory on 
2003 June 23-24 in order to probe farther out into the 
LSB disk of the galaxy. NGC 5949 was imaged in B, 
V, R, and / at the Nickel on 2003 September 11, and 
NGC 6689 was observed in Sloan r' and i' bands with 
the Mosaic camera on the 4 m Mayall Telescope at Kitt 
Peak on 2003 October 20. All imaging took place un- 
der photometric conditions. Reduction and photometric 
calibration of these images followed the description of Si- 
mon et al. (2003). To extend our set of images to the 
near-infrared, we used the 2MASS J, H, and Ks Atlas 
images of each galaxy. Three-color optical images of all 
four galaxies are displayed in Figure 1. 

2.2. Isophotal fits and Stellar Disk Rotation Curves 

We used the optical/near-IR images of these galaxies 
to place reasonable limits on the contributions of their 
stellar disks to their rotation curves. We extracted sur- 
face brightness profiles from the images with the IRAF 
task ELLIPSE, as described in Simon et al. (2003). For 
each image of each galaxy, we ran ellipse allowing the 
center position, position angle (PA), and ellipticity of 
the isophotes to vary with radius. We then defined the 
isophotal center of the galaxy in that band to be the 
weighted average of the ellipse centers. Ellipse was 
then run again with the center fixed, and we measured 
the weighted averages of the PA and ellipticity. We then 



averaged the isophotal centers'"', position angles, and el- 
lipticities for each galaxy in every available band to define 
fiducial values. In cases where a trend as a function of 
wavelength was noticed for any of the parameters, values 
for the redder bands were given preference to minimize 
the importance of extinction. Finally, we ran ELLIPSE 
one more time with all three parameters fixed in order 
to determine the final surface brightness profiles, which 
are shown in Figure 2. The measured isophotal param- 
eters for each galaxy are listed in Table 1. Discussion 
of the various surface brightness profiles can be found in 
§3.5. 

Bell & de Jong (2001) showed that the color and stel- 
lar mass-to-light ratio (M*/L) of a galaxy are correlated. 
In order to obtain an estimate at the mass-to-light ra- 
tios of our targets, we used the observed colors and the 
updated color-M/L relations given by Bell et al. (2003) 
to compute expected R band mass-to-light ratios. For 
NGC 6689, since we do not have any optical observa- 
tions in the Johnson- Cousins system, we instead used 
the relations Bell et al. (2003) defined for the Sloan mag- 
nitude system. Measiired colors and mass-to-light ratios 
are given in Table 2. Two apparent flaws in the Bell et 
al. (2003) models are evident from our calculations: 1) 
the B — V and B — R colors do not predict consistent R 
band mass-to- light ratios, and 2) the predicted K band 
mass-to-light ratios appear to be significantly too high. 

Stellar rotation curves were calculated via numerical 
integration of the stellar surface mass densities (derived 
from the observed surface brightness profiles) using the 
NEMO software package (Teuben 1995). The primary 
assumption required by the NEMO implementation of 
this method is that the disks are infinitesimally thin. 
Removing this assumption leaves the shape of the stellar 
rotation curve unchanged but modestly lowers its ampli- 
tude. 

2.3. Velocity Field Fitting 

The velocity fields of the four galaxies are shown in 
Figure 3. The circles represent fiber-based Ha measure- 
ments, and the fiUed-in regions show CO measurements. 
In general, the velocity fields appear quite regular, with 
only NGC 4605 showing any sign that its kinematic and 
photometric minor axes might not be aligned. 

2.3.1. Comparison of Ha and CO Data 

For NGC 4605 and NGC 5963, we have resolved CO 

velocity fields in addition to the Ha observations. To 
check for agreement between these different data sets, 
we compared the Ha and CO velocities at every posi- 
tion where emission is visible from both species. In the 
case of NGC 5963, we find a mean offset of less than 1 
km s^^, with a scatter of 7.8 km s~^, indicating excel- 
lent agreement between the Ha and CO velocities. For 
NGC 4605, only the westernmost molecular cloud has a 
velocity that agrees exactly with the velocity of the sur- 
rounding ionized gas. Elsewhere in the galaxy there is 
a small but consistent offset of 4.8 km s^^, in the sense 
that the Ha velocities are smaller. The scatter between 

^ For the two galaxies with visually obvious nuclei (NGC 4605 
and NGC 5949), we found that the isophotal centers were located 
quite close to the nuclei. In these cases, we fixed the ellipse centers 
on the nuclei instead of the formal isophotal centers. 
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Fig. 1. — Three-color images of our four target galaxies. Each image uses the asinh scaling recommended by Lupton et al. (2004) to 
retain an optimal combination of color and intensity information, (a) BVR image of NGC 4605 from Lowell Observatory. The southeast 
^ of the galaxy is clearly brighter than the remainder, (b) BVR image of NGC 5949 from Lick Observatory. The blacked-out rectangle 
on the east side of the galaxy is where a large dust grain was not completely removed by fiatfielding. We masked out the affected area for 
the photometric analysis, (c) BVR image of NGC 5963 from Lowell Observatory, (d) r'Hai' image of NGC 6689 from Kitt Peak National 
Observatory. 



TABLE 1 

Measured Isophotal Parameters of Target Galaxies 



Galaxy a (J2000.0) 5 (J2000.0) PA ellipticity inclination 



NGC 4605 


12''39™59!27 


61°36'33'.'3 


119° 


0.630 


71.5° 


NGC 5949 


15'"28™00!69 


64°45'47^'7 


144° 


0.535 


64.6° 


NGC 5963 


15''33'"2r:90 


56°33'35''0 


54° 


0.320 


48.4° 


NGC 6689 


18''34'"50?18 


70°31'27'.'l 


170° 


0.690 


76.0° 



the two velocity fields is 8.4 km s~^. Nevertheless, tilted- 
ring models based only on the CO velocity field match 
Ha-only tilted-ring models quite closely, so this velocity 
offset does not affect our results. Velocity differences of 
this magnitude are expected to arise from flows of the 



ionized gas away from molecular clouds (Fich, Treffers, 
& Bhtz 1982; Fich, Dahl, & Treffers 1990). Since the 
Ha and CO data appear to be both reliable and consis- 
tent with each other, the velocity field fitting described 
in the following section was performed on the combined 
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Fig. 2. — (a) Optical and near-infrared surface brightness profiles of NGC 4605. In each band, the nucleus and exponential disk are 
visible. The black solid line represents an exponential fit to the 7-band profile; the maximum deviation from the data is 0.13 magnitudes, 
(b) Optical and near-infrared surface brightness profiles of NGC 5949. In each band, the nucleus, exponential inner disk, and exponential 
outer disk are all visible, (c) Optical and near-infrared surface brightness profiles of NGC 5963. (d) Optical and near-infrared surface 
brightness profiles of NGC 6689. For all galaxies, the J and Ks data can be traced further out, but we do not plot the data beyond where 
the uncertainties reach a factor of two (0.75 mag). The //-band profiles of each galaxy have been omitted for clarity; the error bars for H 
and Ks overlap at most radii. The vertical dashed lines in panels b), c), and d) emphasize breakpoints in the profiles. The galaxies are 
generally well described by exponential disks, except for NGC 5963, which has an unusual and difficult to characterize surface brightness 
profile. We speculate that this galaxy may have evolved a dense, bulgclikc central region through secular processes (a pseudobulgc). 



velocity data for these two galaxies. 

2.3.2. Tilted-Ring Modeling and Error Analysis 

We derived rotation curves for each galaxy in the same 
way as in Simon et al. (2003). Using the geometric pa- 
rameters measured from the photometry (and listed in 
Table 1), we ran the tilted- ring modeling routine ring- 
fit to extract the rotational, radial, and systemic veloc- 
ities as a function of radius from the velocity fields. As 
detailed in Simon et al. (2003), RINGFIT is similar to 
the standard rotcur algorithm (Begeman 1987) in that 
it includes both rotation and radial motions in its ring 
fits, but it does not allow the PA, inclination angle, or 
center to vary from ring to ring. Next, we applied a boot- 
strap technique to estimate the uncertainties in the PA, 
center position, and inclination angle. Typical bootstrap 
uncertainties were 1" in RA and Dec, 4° in PA, and 2° 
in inclination. For all four galaxies, we found that the 
kinematic and photometric values of the PA, center, and 
inclination angle agreed within the errors. 

After running the bootstrap, we constructed Gaussian 
distributions of each of the geometric parameters, assum- 
ing that the full width at half maximum of the Gaussian 
was equal to the bootstrap uncertainty. We then created 



1000 Monte Carlo realizations of each rotation curve by 
randomly selecting a PA, center position, and inclina- 
tion angle from the Gaussian distributions and running 
RINGFIT with the selected parameters as inputs. This 
process resulted in 1000 Monte Carlo rotation velocities 
at each radius. We then defined our final rotation curves 
to be the means of the Monte Carlo rotation velocities, 
and the uncertainties on the rotation curves to be the 
standard deviations of the Monte Carlo rotation veloci- 
ties. The final radial velocity curves were derived in the 
same manner. By incorporating the uncertainties of the 
PA, inclination, and center position, this technique yields 
more realistic rotation curve error bars than simply prop- 
agating the very small velocity uncertainties from each 
spectrum through the analysis. The Monte Carlo uncer- 
tainties (which we refer to as systematic uncertainties) 
are always much larger than the statistical uncertainties 
calculated by standard error propagation. The uncer- 
tainties we use for the remainder of the paper are the 
sum in quadrature of the systematic and statistical un- 
certainties. 



2.3.3. Higher Order Harmonic Fits 
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TABLE 2 
Stellar Mass-to-Light Ratios 



Mean Inner M,/Lii M*/Lx 
Disk Color^ [Mq/Lq.a] [Mq/Lq^k] 



NGC 4605 

Predictions from B — V color'' 
Predictions from B — R color'' 
Observed maximum disk values 

NGC 5949 

Predictions from B — V color'' 
Predictions from B — R color'' 
Observed maximum disk values 

NGC 5963 

Predictions from B — V color'' 
Predictions from B — R color'' 
Observed maximum disk values 

NGC 6689 

Predictions from r' — i' color'' 
Observed maximum disk values 



0.45'' 0.94 0.72 

0.82" 1.09 0.71 

1.36 0.37 

0. 63'' 1.48 0.76 

1. W^ 1.80 0.78 

2.07 0.49 

0.51" 1.09 0.73 

0.97" 1.38 0.74 

2.09 0.61 

0.22f 1.968 0.78 

3.14 1.07 



''^These colors have been corrected for Galactic extinction and internal ex- 
tinction. The Galactic extinction corrections were taken from Schlegel et al. 
(1998) and the internal extinction corrections were derived with the method 
of Sakai et al. (2000) using the axis ratios from our isophotal fits and H I 

velocity widths extracted from LEDA. 
''Calculated from the relations given in Appendix A of Bell et al. (2003). 
=NGC 4605 colors were measured for 10" < r < 140". 
''NGC 5949 colors were measured for 6" < r < 32". 
"NGC 5963 colors were measured for 8" < r < 32". 
'NGC 6689 colors were measured for 16" <r < 138". 

^Note that this mass-to-light ratio is actually calculated for the Sloan r' 
band, not the Kron-Cousins R band. 



In addition to the standard use of ringfit described 
in the previous section, the algorithm can also be used 
to decompose the velocity field into Fourier components 
up to third order (sin 30 and cosW). If present, these 
higher order terms strongly suggest the existence of a 
bar or other non-axisymmetric perturbation to the grav- 
itational potential. With enough S/N, the unique signa- 
tures of a bar, spiral arms, an elliptical potential, or lop- 
sidedness can be detected (Schoenmakers, Franx, & de 
Zeeuw 1997). We find weak evidence for nonzero higher 
order components in the velocity fields of NGC 4605, 
NGC 5949, and NGC 5963. Each galaxy has several adja- 
cent rings in which the fits deviate from zero by ^ 2a (all 
four higher-order terms for NGC 4605, cos2^^ and sin 30 
for NGC 5949, and cos 261 and sin 29 for NGC 5963). Typ- 
ical amplitudes of the higher order components are 2 — 5 
km s~^, which should not be large enough to affect our 
derivation of the density profiles in the following section. 
Although the significance of these components docs not 
generally exceed 2.5c7 in a given ring, we argue that these 
detections are probably real because 1) a number of con- 
secutive rings have consistent Fourier components, and 
2) the rotation curve uncertainties appear to be overesti- 
mated (see §3.3). NGC 6689 docs not have enough data 
points in most of its rings to obtain reliable higher-order 
fits. For comparison, we find 3 — 5 km s~^ cos 29, cos 30, 
and sin 30 components over small ranges of radii in the 



NGC 2976 velocity field studied by Simon et al. (2003). 
Including terms up to third order in the tilted-ring mod- 
els (seven free ringfit parameters instead of three) does 
significantly improve the reduced vahie of the fits in 
many cases, but this is often because the number of data 
points in the ring is not much larger than seven so that 
the fit can go through every point. More detailed de- 
scriptions of the higher order motions in each galaxy can 
be found in §3.5. 

3. RESULTS 

In this section we fit the rotation curves with a variety 
of functional forms to ascertain the shapes of the density 
profiles. We first consider the likely contribution of the 
stellar disk to the rotation curve, and then remove it from 
the data to reveal the rotation curve of the dark matter 
halo. 

The rotation curves produced by the tilted-ring models 
in §2.3.2 are displayed in Figures 4a - 7a. Radial velocity 
curves and the uncertainties on both the rotation veloci- 
ties and the radial velocities are also plotted in the same 
figures. 

3.1. Removing the Baryons 

What we are actually interested in are the rotation 
curves of the dark matter halos of these galaxies, which 
are only equivalent to the total rotation curves shown in 
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Fig. 3. — Velocity fields of the four galaxies discussed in this paper, (a) Combined Ha (circles) and CO (fiUed-in) velocity field of 
NGC 4605. (b) Ho velocity field of NGC 5949. (c) Combined Ha and CO velocity field of NGC 5963. (d) Ho velocity field of NGC 6689. 



Figures 4a - 7a in the minimum disk scenario (Af,/L = 
0). We therefore need to account for the contribution of 
the baryons to the rotation curves. An upper limit on the 
dark matter rotation curve (and also the slope of the den- 
sity profile) can be found if the disk mass is zero, and a 
lower limit to the dark matter rotation curve and density 
profile slope is obtained for a maximum disk. In general, 
for galaxies of normal surface brightness, the minimum 
disk solution is physically unrealistic and the actual mass 
distribution is likely to be closer to the maximum disk 
case. 

We used the stellar disk rotation curves calculated 
in §2.2 to construct maximum disk models^. In order 
to avoid being unduly influenced by a single unrepre- 
sentative point in any of the observed rotation curves, 
we fit the inner five points of each rotation curve with 
a smooth function (a power law). We then incremen- 
tally increased the stellar mass-to-light ratio (the stel- 
lar rotation curve scales ^jK-h/L) until the 
stellar rotation curve exceeded the fit to the observed 

With the possible exception of NGC 6689, these galaxies do 
not appear to contain bulges, so the stellar disk is the only baryonic 
component that we model. Including the gas disk as well would 
increase the effective mass-to-light ratio of the disk somewhat, but 
does not affect our results (see following paragraph). 



points. Because of the influence of the nucleus on the 
stellar rotation curve of NGC 4605, we ignored the in- 
nermost point of the fit to the observed rotation curve 
for that galaxy. With this procedure we defined the fol- 
lowing maximum disk mass-to-light ratios: M^,/L]^ = 
IMMq/Lqr (NGC 4605), mJLr = 2.07Mq/Lqr 
(NGC 5949)', M^/Lr = 2.O9Mo/L0 r (NGC 5963), and 
AU/Lr' = S.WMq/Lq,^/ (NGC 6689). Typical uncer- 
tainties on these values are ^ 15 %. For all of the galax- 
ies, these kinematic mass-to-light ratios are significantly 
higher than the mass-to-light ratios predicted from the 
galaxy colors by the Bell et al. (2003) population synthe- 
sis models (see Table 2). We conclude the disks of these 
galaxies are probably somewhat submaximal. For the 
remainder of this paper, we adopt the following mass-to- 
light ratios (the average predicted from the B — V and 
B-i? colors shown in Table 2): M^/Lr = 1MMq/Lq r 
(NGC 4605), M,/Lr = IMMq/Lqr (NGC 5949), 
AU/Lr = 1.24:Mq/Lqr (NGC 5963),' and M^/L^- = 
1.96Mq/L0,^. (NGC 6689). 

Note that since we do not have H I data, we are ne- 
glecting the contribution of the gas to the rotation curve 
in this analysis. In galaxies the size of our targets, how- 
ever, the stellar disk almost always contributes signifi- 
cantly more mass at the center of the galaxy than the 
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Fig. 4. — (a) Tilted-ring model for NGC 4605. The thick blue curve represents the rotation curve, and the thick red and orange curves 
show the radial velocity curve and the systemic velocities, respectively. Starting at a radius of 72" the systemic velocities were held fixed 
at zero. The shaded gray regions surrounding the rotation and radial velocity curves represent the combined systematic and statistical Icr 
uncertainties from the bootstrap and Monte Carlo analysis. Note that for convenience we have plotted the radial velocities as positive; 
this choice does not indicate whether the radial motions are directed inwards or outwards. The thin black curve is a power law fit to the 
rotation curve. The residuals from this fit are plotted in the lower panel, and the Icr and 2a scatter of the data points around the fit is 
shown by the shaded gray areas, (b) Disk-subtracted rotation curve of NGC 4605. The gray points are the original rotation velocities from 
panel (a), and the black points are the dark matter halo rotation curve after subtracting the stellar disk (M,/Lfl = LOIMq/Lq^a), which 
is shown as a gray dashed line. The thick green, cyan, and magenta curves show power law, NFW, and pseudo-isothermal fits to the halo 
rotation curve, respectively. The residuals from these fits are displayed in the lower panel, and the \a and 2a scatter of the data points 
around the power law fit is shown by the shaded gray areas. 




gas does and the effect of including the gas is similar to 
a 20% change in the stellar mass-to-light ratio (Bolatto 
et al. 2002; Simon et al. 2003). If we arbitrarily in- 
crease M*/L to simulate this effect, the slope of the dark 



matter density profile (see below) decreases by 2 — 12%. 
Allowing the stellar disk to have a nonzero thickness (see 
§2.2) would offset this decrease. 
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3.2. Power Law Density Profile Fits 

We can now consider the central density profile slopes 
of the dark matter halos of each galaxy. Following the 
notation of Simon et al. (2003), we parameterize the 
rotation curve by Vrot oc 7-(^^"°")/^, which corresponds 
to a power law density profile: p^M oc r~"°". Allow- 
ing the stellar mass-to-light ratios to vary from zero up 
to the maximum disk values determined in the previous 
paragraph, the limits we place on aum are given in Table 
3. Obviously, the impact of changing the mass-to- light 
ratio on the density profiles is minimal for NGC 4605 



and NGC 5949. These galaxies have very similar dark 
matter density profiles in the maximum and minimum 
disk cases, so their density profile slopes are tightly con- 
strained. 

The rotation curve fits for these limiting cases are dis- 
played in Figures 4-7. In the left panel of each figure, a 
minimum disk power law fit is displayed, with the resid- 
uals from the fit in the lower panel. The right panel of 
each figure contains a power law fit to the dark mat- 
ter halo in the maximum disk case, again with residuals 
shown in the lower panel. In all cases a power law pro- 
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TABLE 3 

Limits on Dark Matter Density Profile 

Slopes 



Galaxy maximum disk minimum disk 



NGC 2976 


0.01 


± 


0.13 


0.27 


± 


0.09 


NGC 4605 


0.71 


± 


0.06 


0.90 




0.02 


NGC 5949 


0.79 


± 


0.17 


0.93 


± 


0.04 


NGC 5963 


0.75 


± 


0.10 


1.41 




0.03 


NGC 6689 


0.43 


± 


0.18 


1.07 


± 


0.06 



vides a reasonable fit to the rotation curve at least out to 
the radius at which the rotation curve begins to flatten 

significantly (at which point the density profile slope is 
obviously changing with radius and a single power law 
cannot be expected to accurately represent the rotation 
curve). The numerical parameters of the power law fits 
for the disk-subtracted rotation curves are listed in Ta- 
ble 4, and for the total rotation curves (with no disk 
subtraction) in Table 5. The dark matter density profile 
slopes we derive after subtracting the stellar disks are 
aDM = [0.78, 0.88, 1.20, 0.79] for NGC 4605, NGC 5949, 
NGC 5963, and NGC 6689, respectively. If we fit the 
velocity fields without allowing radial motions, we find 
density profile slopes of a^M = [0.64, 0.98, 1.19, 0.79]. 
These changes are small (comparable to or less than the 
1(7 uncertainties) and go in different directions for differ- 
ent galaxies. We conclude that including radial motions 
in our tilted-ring models does not systematically affect 
the density profile slopes we derive. 

3.3. Alternative Density Profile Fits 

In order to obtain a more complete understanding of 
the allowed shapes of the inner density profiles of the 
galaxies in our sample, we also fit the rotation curves 
with several other functional forms. The two most com- 
monly used profiles in the literature are the pseudo- 
isothermal profile (an isothermal sphere with a central 
core) and the NFW profile^. The functional form for the 
pseudo-isothermal density profile is 

where po is the central density and Vc is the core radius 
(e.g., Kent 1986). 

Since pseudo-isothermal density profiles are not ex- 
pected on theoretical grounds, and there is little obser- 
vational support for NFW profiles, it is important not 
to limit ourselves to only these possibilities. Fitting the 
data with other density profiles can help to elucidate the 
physical meaning of the fits. We therefore also used a 
profile we defined to have characteristics intermediate be- 
tween those of the pseudo-isothermal and NFW profiles: 

^("■^ = (r/re)V2(l + r/re)3/2 ' 

where rc and pc are the characteristic radius and den- 
sity. The inner slope of this profile is a = 0.5, giving it a 

For our NFW fits, we chose the scale radius Vs and the virial 
radius r2oo to be the two free parameters. 



shallow inner cusp. At large radii the slope is q = 2 (cor- 
responding to a fiat rotation curve), so the fact that our 
observations only just reach the flat part of the rotation 
curve should not adversely affect the fits. Note that the 
values for fits to this function are given in Tables 4 - 
5 even though we do not list rc and pc for each galaxy. 

We used a ^^-minimization routine to find the best 
fit parameters for each of these functions. With the ex- 
ception of NGC 4605, in which small-scale bumps and 
wiggles in the rotation curve inflate x^) typical reduced 
values for the best fits were w 0.3, suggesting that 
the systematic uncertainties on our rotation curves have 
been overestimated by a factor of about 1.7. Button et 
al. (2003) similarly inferred that our rotation curve er- 
ror bars for NGC 2976 were overestimated, indicating 
that the bootstrap and Monte Carlo analysis described 
in §2.3.2 may be excessively conservative. 

The results of these fits after removing the contribu- 
tion of the stellar disk of each galaxy are summarized in 
Table 4, and fits without removing the stellar disk ro- 
tation velocities are given in Table 5. The best power 
law, NFW, and pseudo-isothermal fits (after subtracting 
the stellar disk) are displayed in Figures 4b - 7b by the 
green, cyan, and magenta curves, respectively. Residuals 
from each fit are shown in the same colors in the lower 
right panels of each figure. 

With the exception of NGC 5963 (and NGC 6689 in 
the minimum disk case), we were unable to obtain rea- 
sonable NFW fits to the rotation curves. Since our data 
generally do not extend well into the flat part of the 
rotation curve, Vg and r2oo are highly covariant. To pre- 
vent both parameters from running away to unphysically 
large values for NGC 4605, NGC 5949,' and NGC 6689, 
we held fixed and only fit for r2oo • An infinite number 
of such fits are possible, with decreasing as and r2oo 
increase. The NFW fits listed in Tables 4 and 5 are for 
rs values chosen to keep x^ low (Ax^ ~ 1 relative to the 
best fitting power law) without allowing r2oo to become 
completely unreasonable. 

3.3.1. Fit Results 

NGC 5963, strikingly, is best fit by a very steep (q;dm > 
1) power law or an NFW profile unless its stellar mass- 
to-light ratio is unrealistically high. The formal best fit 
for NGC 4605 is a pseudo-isothermal profile, although 
none of the profiles have good x^ values. The best power 
law (which does not describe the full rotation curve very 
accurately) still has a relatively steep slope. NGC 5949 
is also best fit by pseudo-isothermal density profile, but 
power laws with slopes slightly shallower than NFW, and 
the modified pseudo-isothermal profile given in Equa- 
tion 2, fit very well. NGC 6689 is better described by 
a power law than a pseudo-isothermal profile, and in the 
case of a low stellar M/L can even be; fit by an NFW 
profile, although the scale radius and virial velocity sug- 
gested by the NFW fit are unrealistically large. For ev- 
ery galaxy except NGC 4605, Equation 2 (the modified 
pseudo-isothermal profile with a shallow central cusp) 
provides fits that are comparable to or better than the 
pseudo-isothermal profile. 

3.3.2. Navarro et al. (2004) Density Profile 

In addition to the profiles described above, we also fit 
the rotation curves (after removing the stellar disks) with 
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TABLE 4 

Rotation Curve Fit Results After Disk Subtraction 



Power law NFW Pseudo-isothermal 

fit parameters fit parameters fit parameters 



Galaxy 


Vb 


a 


re 


^200 


c 


V200 


PO 




Rc 




x'nfw'^ 


2 d 


w ■ 




[km s-i] 




[kpc] 


[kpc] 




[km s-l] 


[g cm- 




[kpc] 








NGC 2976* 


43.0 


0.01 


... b 


... b 


... b 


... b 


7.7 X 10 


-24 


2.5 


0.27 


> 5.8 


0.43 


1.00 


NGC 4605 


51.4 


0.78 


> 5.00 


> 91 


< 18.3 


> 65 


2.0 X 10 


-23 


0.8 


3.74 


< 8.03 


1.53 


2.88 


NGC 5949 


39.3 


0.88 


> 6.40 


> 87 


< 13.6 


> 62 


7.8 X 10 


-24 


1.2 


0.38 


< 1.24 


0.24 


0.27 


NGC 5963 


66.4 


1.20 


10.8 


160.9 


14.9 


114.3 


3.4 X 10 


-23 


0.8 


0.34 


0.43 


1.36 


0.64 


NGC 5963"= 


74.5 


1.50 


6.71 


127.7 


19.0 


90.7 


2.7 X 10 


-23 


0.9 


2.19 


0.42 


0.98 


0.57 


NGC 6689 


37.8 


0.79 


> 8.7 


> 108.8 


< 12.5 


> 77 


4.8 X 10 


-24 


2.1 


0.30 


< 1.31 


0.46 


0.30 



^Using data from Simon ct al. (2003). 

''No meaningful NFW fits could be made for NGC 2976, even with only one parameter free. The "best" fit had Xred ~ 

•^Using the rotation curve presented in this paper plus the rotation velocities measured by Bosma et al. (1988) at 55" < r < 165". 

''All values listed here actually represent the reduced (x^ per degree of freedom in the fits). 



TABLE 5 

Rotation Curve Fit Results With No Disk Subtraction 



Power law NFW Pseudo-isothermal 

fit parameters fit parameters fit parameters 



Galaxy 


Vb 
[km s-i] 


a 


[kpc] 


^200 

[kpc] 


c 


V200 
[km s-l] 


PO 

[g cm-3] 


Rc 
[kpc] 


2 d 
Xpt 


2 d 
Xnfw 


2 d 




NGC 2976* 


52.1 


0.27 


... b 


... b 


... b 


... b 


1.3 X 10-23 


1.5 


0.33 


> 5.8 


0.48 


0.93 


NGC 4605 


65.0 


0.89 


> 42.0 


> 347 


< 8.3 


> 246 


4.0 X 10-23 


0.7 


8.28 


< 9.28 


1.53 


4.92 


NGC 5949 


63.3 


0.93 


> 38.0 


> 325 


< 8.6 


> 231 


2.2 X 10-23 


1.1 


1.50 


< 2.00 


1.00 


0.64 


NGC 5963 


87.8 


1.41 


2.86 


102.7 


35.9 


72.9 


1.3 X 10-22 


0.4 


6.88 


3.57 


1.22 


3.04 


NGC 5963^^ 


89.5 


1.62 


3.42 


109.9 


32.1 


78.0 


1.2 X 10-22 


0.4 


6.52 


2.64 


1.07 


1.82 


NGC 6689 


53.1 


1.07 


41.5 


286.7 


6.9 


203.6 


1.4 X 10-23 


1.1 


0.32 


0.31 


0.95 


0.42 



^Using data from Simon ct al. (2003). 

''No meaningful NFW fits could be made for NGC 2976, even with only one parameter free. The "best" fit had Xred ~ ^'^^ 
•^Using the rotation curve presented in this paper plus the rotation velocities measured by Bosma et al. (1988) at 55" < r < 165". 
''All values listed here actually represent the reduced x2 (x^ per degree of freedom in the fits). 



the new density profile proposed by Navarro et al. (2004): 



Priir) = P-ie 



(3) 



where p-i and r_2 are the density and radius at which 

the logarithmic density profile slope equals —2 and t] is 
an additional free parameter® that controls how fast the 
density profile slope changes with radius. The rotation 
curve associated with this density profile is 



3 2 



R 



.(4) 

where 7(0,2;) is the lower incomplete gamma function. 
With the extra free parameter afforded by this func- 
tion, we are able to achieve very good fits to the rotation 
curves of each of our five targets. In Table 6 we present 
the results of these fits and compare them to the fits 
Navarro et al. (2004) performed to four simulated dwarf 

* Navarro ct al. (2004) use a instead of r;, but since wc have 
defined a to be the slope of a power-law density profile we adopt 
a diflferent symbol here for clarity. 



galaxies that have total masses similar to our targets. 
Only for NGC 5963 are the fit parameters we derive re- 
motely close to those measured by Navarro et al. (2004) 
in their simulations. 

3.4. Comparison to Previously Published Rotation 

Curve Data 

In the remainder of this section, we discuss details of 
i/gur analysis of each galaxy. Readers who are primarily 
interested in the more general results of our work may 
wish to skip to §4. 

Each of the four galaxies we observed have existing 
Ha or H I rotation curve data in the literature. Our new 
data should improve on previous measurements in veloc- 
ity resolution, control of systematics, angular resolution 
(in the case of H I measurements) , and sensitivity, but it 
is useful as a first-order test of the accuracy of our data 
to compare the rotation curves we obtain with previous 
measurements. 

NGC 4605 has been the subject of several rotation 
curve studies (Rubin, Thonnard, & Ford 1980; Sofue et 
al. 1998; Bolatto et al. 2002). All three of those stud- 
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TABLE 6 

Navarro et al. (2004) Profile Fit Results 



Galaxy 


r-2 




P 


-2 




V 


^2a 




[kpc] 


[g 


cm ' 








NGC 2976'' 


1.65 


5.6 


X 


IQ- 


-24 


8.99 


0.26 


NGC 4605 


1.34 


5.8 


X 


IQ- 


-24 


1.84 


1.37 


NGC 5949 


3.85 


7.5 


X 


10" 


-25 


0.67 


0.25 


NGC 5963 


6.44 


6.0 


X 


10" 


-25 


0.28 


0.62= 


NGC 5963'^ 


6.04 


6.4 


X 


10" 


-25 


0.28 


0.43 


NGC 6689 


3325 


1.3 


X 


10" 


-28 


0.12 


0.31 





4.55 


1.1 X IQ- 


-25 


0.164 


02"^ 


4.28 


1.5 X IQ- 


-25 


0.211 


03"^ 


3.62 


1.5 X 10- 


-25 


0.122 


04"^ 


3.62 


2.1 X 10" 


-25 


0.166 



'^X^ values listed here actually represent the re- 
duced (x^ per degree of freedom in the fits) 

''Using data from Simon et al. (2003). 

■^The three-parameter fit using only our data for 
NGC 5963 is degenerate. Better values than 
shown here can be achieved, but only when r_2 and 
p—2 run away to uiu'calistic values. We therefore 
fixed r) to the value obtained when the Bosma et al. 
(1988) data are included as well (see the following 
line of this table) and ran the fit with only two free 
parameters. 

Using the rotation curve presented in this paper 
plus the rotation velocities measured by Bosma et 
al. (1988) at 55" <r < 165". 

"Simulated dwarf galaxy halos from Navarro et al. 
(2004). 



ies derived very similar rotation curve shapes and am- 
plitudes. While our methodology is different from that 
employed by the previous authors (we do not produce 
separate fits for the approaching and receding sides of 
the galaxy), our results are qualitatively in agreement. 
In particular, the break in the rotation curve at a radius 
of 40" is quite apparent in the two recent data sets (the 
Rubin et al. 1980 rotation curve docs not have enough 
angular resolution) . and the maximum rotation velocities 
match well. The primary quantitative discrepancy is the 
result of the incorrect plate scale assumed by Bolatto et 
al. (2002) for the Ha spectrum they obtained. In that 
paper, we used a value of 3.0 arcsec pixel^^, but subse- 
quent measurements showed that the actual plate scale 
is 2.1 arcsec pixel" ^. Correcting this error reduces the 
extent of the Ha rotation curve by 30%, but should not 
substantially change the power law indices measured for 
the rotation curve and density profile. Indeed, our new 
rotation curve matches the old one near the center of the 
galaxy (r < 10") and at large radii (r > 65"). From 
10" to 30" our rotation curve is a few km lower than 
the old one, and from 30" to 65" our rotation curve ex- 
ceeds the old one by 1 — 6 km s~^. The net effect of 
these changes^ is to modestly steepen the maximum disk 
density profile from Pdm oc r 

(Bolatto et al. 2002) 

® Note that the stellar disk we use is not the same as that of 
Bolatto et al. (2002). Bolatto et al. calculated the stellar 
rotation curve for an exponential disk and we use the actual surface 
brightness profile to determine the stellar rotation curve. Although 
the resulting stellar disks are not quite identical, the differences do 
not significantly affect the density profile results. 



to Pdm r . Given typical uncertainties on aoM of 
0.1 in our analysis (Simon et al. 2003), these results 
are marginally consistent. The earlier study by Persic 
& Salucci (1990) found a density profile of p oc r~°'^^, 
also consistent with our results. The best fitting power 
law rotation curve at large radii still overestimates the 
rotation velocities in the inner region of the galaxy, as 
seen in Bolatto et al. (2002) and Figure 4a. 

The rotation curve of NGC 5949 has previously been 
measured by Karachentsev & Petit (1990) and Courteau 
(1997) with long-slit spectra. The Karachentsev & Petit 
data are consistent with the rotation curve we derive, 
although our rotation velocities are somewhat larger and 
our measurements extend to larger radii. Our rotation 
curve matches that of Courteau within the uncertainties 
of the data. 

The kinematics of NGC 5963 have been investigated 
in detail in two previous papers (Romanishin, Strom, 
& Strom 1982; Bosma, Athanassoula, & van der Hulst 
1988). The shape of the Ha rotation curve measured by 
Romanishin et al. (1982) is very similar to the shape of 
our rotation curve; the amplitude they derive is slightly 
larger (asymptotic rotation velocity of 131 km s~^), but 
this appears to be the result of the smaller inclination 
angle they used. The H I rotation curve presented by 
Bosma et al. (1988) is perfectly consistent with our rota- 
tion curve over the common range of radii, except that 
the inner two points of their rotation curve are affected 
by beam smearing. For the fits presented in Tables 4 - 
6, we give results for both our data alone, and our data 
with the Bosma, Athanassoula, & van der Hulst (1988) 
points at large radii added. 

The Ha velocity field of NGC 6689 was observed by 
the GHASP survey (Garrido et al. 2003). Those authors 
found a decidedly asymmetric rotation curve, but the 
filter they used for the observations may have cut off 
the emission line profiles on the approaching side of the 
galaxy. The receding side of the rotation curve presented 
by Garrido et al. (2003) is consistent with the overall 
rotation curve we derive out to a radius of 45", but at 
larger radii our rotation curve shows significantly smaller 
rotation velocities. 

3.5. Comments on Individual Galaxies 

3.5.1. NGC 4605 

Of our four target galaxies, NGC 4605 is the only one 
that appears to deviate significantly from axisymmetry. 
The non-axisymmetric structures are visible both photo- 
metrically and kinematically. In broadband optical im- 
ages, the galaxy contains a pronounced elongated region 
that is offset from its nucleus to the east. This feature, 
which has a somewhat higher ellipticity and a different 
position angle than the rest of the galaxy, persists out to 
K-band although it weakens with increasing wavelength. 
In Figure 8, we plot the isophotal fit parameters as a 
functicm of radius (where we now allow the isophotal cen- 
ter, PA, and ellipticity to vary from ring to ring), and the 
impact of the asymmetry is easily visible out to a radius 
of 70". Since this lopsidedness is also apparent in our 
Ha image (the galaxy's Ha emission extends to twice 
as large a radius on its eastern side as on the western 
side) and CO map, it may be associated with recent star 
formation. 
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Fig. 8. — il-band isophotal fit parameters for NGC 4605 with the eUipticity, PA, and center left free to vary, (a) Isophotal eUipticity 
as a function of radius. The peaJs at r = 30 — 40" is due to the asymmetric structure of the galaxy. Note the corresponding features in 
panels (b) and (c). The horizontal line shows the mean value used in our analysis, (b) Isophotal position angle as a function of radius. 
The horizontal line shows the mean value used in our analysis, (c) Isophotal center position along the x (right ascension) axis of the image, 
in arcseconds. The position of the nucleus corresponds to Xq = 0, which coincides with the center of the outermost isophotes. The center 
position along the y (declination) axis shows almost identical behavior, but with smaller amplitude because the asymmetric structure is 
located mostly east of the galaxy center, (d) Surface brightness as a function of radius, using the parameters shown in the other panels. 
Note that since our final isophotal fits did not allow the parameters to vary with radius, this surface brightness profile differs slightly from 
the one used to derive the stellar disk rotation curve (Figure 2a). 



The kinematic effects of this asymmetry are twofold: 
1) the kinematic center of the galaxy is offset from the 
nucleus, and 2) the galaxy contains relatively large (and 
significant) noncircular motions. Tilted-ring models of 
the galaxy centered on the nucleus result in a system- 
atic trend of the systemic velocity with radius. Moving 
the center position around so as to remove this trend 
and simultaneously minimize the scatter of the systemic 
velocities about their mean value, we found that the pre- 
ferred kinematic center is southeast of the nucleus by 
10". Offsets of the kinematic center of this magnitude 
do not significantly change the slope of the density pro- 
file. Irrespective of the center position, the models also 
show a clear detection of radial motions in NGC 4605. 
The amplitude of those motions increases to 8 km s^^ at 
a radius of 20", and then stays constant out to 80". At 
larger radii (the outermost six rings) the radial motions 
begin to increase again, reaching 20 km s~^ at the edge of 
the galaxy (r = 113"). However, this increase is not very 
well constrained due to the poor angular coverage of the 
velocity data at large radii. Unsurprisingly, we find small 



(2 — 6 km s~^) but consistently nonzero cos 261, sin 2^, 
cos 36, and sin 3^ terms in the velocity field of NGC 4605 
at various ranges of radii between 16" and 68". 

The combination of photometric and kinematic lopsid- 
cdness suggests that the mass distribution of the galaxy 
indeed is lopsided, as opposed to the possibility of asym- 
metric star formation mentioned earlier. Despite the lop- 
sidedness of NGC 4605, the stellar disk is described well 
by an exponential profile, as shown in Figure 2a. 

In order to determine the potential impact of these 
asymmetries on our density profile analysis, we consid- 
ered the maximum changes they could cause. The photo- 
metric structure makes it difficult to choose an eUipticity 
and position angle that are representative of the whole 
galaxy. In §2.2 we used a weighted mean of the eUipticity 
and PA at radii between the outer edge of the asymmet- 
ric structure (at r = 70") and the edge of the galaxy. As 
an alternative, we tried setting the eUipticity to its maxi- 
mum and minimum values over the whole galaxy, and re- 
peated the tilted-ring modeling and isophotal fitting. We 
found that the rotation curve and density profile shape 
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are largely insensitive to changes in the inclination an- 
gle. Even for ellipticitics as small as 0.56 {i = 66°) and 
as large as 0.73 (i = 79°), the largest change in the de- 
rived value of is —0.05. We therefore conclude that 
our analysis is robust to the uncertainties introduced by 
the asymmetry of NGC 4605. 

3.5.2. NGC 5949 

The surface brightness profiles of NGC 5949 are very 
similar to those of NGC 2976 (Simon ct al. 2003). Both 
galaxies contain a nucleus, a shallow (large scale length) 
inner disk, and a steep (small scale length) outer disk. 
Although other disk galaxies arc known to have similar 
structures (e.g., Pohlcn et al. 2002), the origin and physi- 
cal significance of those two apparently distinct disks are 
not understood. Nevertheless, there is no reason to doubt 
that our thin-disk approximation is valid for this galaxy, 
so the shape of the stellar rotation curve we calculated 
should be accurate. Due to the high level of symmetry 
of NGC 5949 in all bands, the photometric inclination 
and PA are extremely well constrained, leading to very 
small systematic uncertainties on the rotation curve. It 
is worth noting that NGC 5949 is the only galaxy in our 
sample in which the shape of the stellar rotation curve is 
very similar to the shape of the observed rotation curve. 
Therefore, it is possible to model this galaxy with essen- 
tially only a stellar disk and barely any contribution from 
dark matter. The stellar mass-to-light ratio required to 
accomplish this is M^/Ln « 2.7Mq/ Lq,ji, 30 % higher 
than the derived maximum disk value. 

Our tilted-ring modeling of NGC 5949 reveals a small 
radial component to the velocity field with modest sig- 
nificance. These radial motions are detected beginning 
at a radius of 20", and remain present out to the edge of 
the galaxy. The maximum amplitude is 8 km s^^, and 
there are a total of six consecutive rings that have ra- 
dial motions deviating from zero by at least la. Since 
the radial motions are < 10% of the rotation speed at 
all radii, they should not affect our mass models of the 
galaxy. Note that because we do not know which side of 
the galaxy is closer to us, we cannot determine whether 
the detected radial motions correspond to inflowing or 
outflowing material. We also detect a --^ 2 km cos 2^ 
term from 24" to 36" and a ^ 3 km sin 3^ term from 
12" to 32". 

3.5.3. NGC 5963 

Unlike the other three galaxies, the luminous compo- 
nent of NGC 5963 does not contain an easily identiflable 
exponential disk. At the center of the galaxy is a bright, 
elongated bar-like feature about 4" across. Outside this 
source is a small disk-like region (500 pc in radius) with 
four tightly wound spiral arms. At a radius of 15" (950 
pc) the surface brightness profile begins a steep decline, 
falling by nearly 3 magnitudes over 18". Surrounding 
this region is an LSB, nearly exponential disk that ex- 
tends out to a radius of at least 120" (see Figure 2c). 
The inner spiral arms can be traced out into this outer 
disk at very low surface brightness levels. 

It is difficult to interpret this surface brightness profile 
in terms of the standard model of a disk galaxy. Ordi- 
narily, one might assume that the bright central region 
of NGC 5963 is a bulge that just happens to be at the 



center of an unusually faint disk. The clear presence of 
spiral arms in this region, however, suggests that oven 
here the galaxy is highly flattened. One possibility is 
that this structure is a pseudobulge that has formed via 
secular evolution of the galaxy (e.g., Kormendy & Ken- 
nicutt 2004) . The most sensible way to derive the stellar 
rotation curve in this case is to apply our standard as- 
sumption that the disk can be treated as infinitesimally 
thin. The assumption that the stellar mass-to-light ra- 
tio remains constant throughout this galaxy may not be 
correct, but since the stars are far from dominating the 
gravitational potential for reasonable mass-to- light ratios 
this should not substantially change our results. 

NGC 5963 also contains noncircular motions. The ra- 
dial term in the tilted-ring fits deviates from zero by more 
than 1 cr in five rings, from 12" to 30". The radial mo- 
tions reach as high as 15 km s^^ in this region, before 
beginning to oscillate about zero (with large uncertain- 
ties) at larger radii. The second-order Fourier terms in 
this galaxy both have amplitudes of ~ 3 km s^^ from 
radu of 16" - 24" (sin 26') and 12" - 24" (cos 26*). 

3.5.4. NGC 6689 

The high inclination of NGC 6689 makes our images 

and surface brightness profiles less revealing of details of 
its structure. Nevertheless, it is clear that this galaxy 
contains an exponential disk and a central light excess 
that could be attributable to a bulge or a bar. No bar 
is evident in our images, but as we noted in §2.1, some 
catalogs classify the galaxy as barred. 

Despite the large inclination angle, the observed Ha 
emission lines from NGC 6689 are well represented by 
Gaussian fits. This implies that sightlines through 
NGC 6689 generally only intersect a single H II region, 
which is consistent with the appearance of the galaxy in a 
narrow band Ha image. Therefore, it is not necessary to 
use an envelope-tracing method (e.g., Sofuc 1996: Gentile 
et al. 2004) instead of Gaussian fits to extract rotation 
velocities from the spectra. 

Neither our standard tilted-ring model nor the higher 
order harmonic fits detected any deviations from circu- 
lar rotation in NGC 6689. To some degree, this is likely 
due to the inclination of the galaxy, which limits our re- 
solving power along the minor axis where radial motions 
are most prominent. Nevertheless, a radial component 
of 10 km s~^ as we found in the other three galaxies 
should have been detected if present. 

4. DISCUSSION 

4.1. Is There a Universal Density Profile? 

The primary goal of this study is to determine as ac- 
curately as possible the inner density profiles of the dark 

matter halos of the galaxies we observed. At radii of < 1 
kpc, does the dark matter density continue to increase, as 
predicted by NFW and numerous other theoretical stud- 
ies? Or are the density profiles flat, with constant den- 
sity cores, as most previous observers have concluded? 
Are the variations in density profile shape from galaxy 
to galaxy as small as the simiilations suggest? 

4.1.1. Observational Results 

For two of the galaxies in our sample (note that for 
the remainder of the paper we include NGC 2976 in our 
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analysis so that wc have a total sample of five galaxies), 
these questions are easy to answer. NGC 2976, as Simon 
et al. (2003) showed, unambiguously contains a constant 
density core. An Udm = power law (constant density) 
provides an excellent fit to the rotation curve. A pseudo- 
isothermal profile with a core comparable to the optical 
size of the galaxy also produces a good fit. NGC 5963, 
on the other hand, quite clearly has a very steep central 
density profile. A power law with a slope of a^^ — 1-20 
fits the rotation curve very well, and an NFW fit with 
rs « 11 kpc and a concentration parameter of 14.9 is 
nearly as good. A pseudo-isothermal fit is significantly 
inferior. 

For the remaining three galaxies, the interpretation of 
the fit results is not nearly as straightforward. NFW 
fits to the disk-subtracted rotation curves can be carried 
out, but the fit parameters are not usefully constrained. 
Power laws with slopes between Q!om = 0.78 (NGC 4605) 
and ttoM = 0.88 (NGC 5949) fit the rotation curves well, 
but they are matched in each case by pseudo-isothermal 
profile fits with comparable or better reduced values. 

How can a rotation curve be simultaneously consis- 
tent with both a pseudo-isothermal halo and an a ~ 0.8 
power law density profile? As shown by van den Bosch 
& Swaters (2001), the difference between the rota- 
tion curves associated with these density profiles tends 
to be smaller than typical observational uncertainties, 
even with the high velocity resolution of our data. The 
pseudo-isothermal profile has a slope that varies con- 
tinuously from (constant density) at its center to -2 
(isothermal) at large radii. A power law, of course, has a 
constant slope that (for our fits) is intermediate between 
these two values. That both profiles fit the data suggests 
that 1) the mean density profile slope over the observed 
region is equal to the value preferred by the power law, 
but that 2) the slope changes with radius, which is bet- 
ter described by the pseudo-isothermal profile. Previous 
studies have often assumed that if a pseudo-isothermal 
rotation curve fits the data, then the galaxy in question 
must contain a constant-density core. Our successful fits 
with the modified pseudo- isothermal profile (Equation 2) 
show that this assumption is not correct. Even when a 
pseudo-isothermal profile fits well, other density profiles 
with steeper central slopes may provide equally good fits. 

Taken together, these observational results argue 
against the proposition that all galaxies share a universal 
density profile. Common two-parameter profiles other 
than a power law cannot fit all of the galaxies in our 
sample, and even for power law fits the solutions span 
the range from constant-density to very cuspy. If a uni- 
versal density profile exists, the scatter from halo to halo 
is large. 

4.1.2. Simulation Results 

The idea of a universal dark matter density profile 
traces back to NFW, and the most recent simulations 
continue to support this picture (Navarro et al. 2004; 
Diemand, Moore, & Stadel 2004: Stoehr 2004). The 
favored functional form of the universal profile, however, 
has not been agreed upon. For several years the de- 
bate focused on the exact value of the central cusp slope, 
but with the increasing resolution of the simulations it 
now seems that this question may be the wrong one to 
ask. The highest resolution ACDM simulations reveal 



that dark matter density profiles do not converge to an 
asymptotic central slope (Power et al. 2003; Navarro et 
al. 2004; Diemand et al. 2004). Instead, the logarithmic 
slope continues to get shallower at smaller radii, and the 
best procedure for extrapolating the profiles to radii be- 
low the resolution limits^" of the simulations is not clear. 

The new density profile proposed by Navarro et al. 
(2004) (Equation 3) fits simulated dark matter halos 
more accurately and over a wider radial range than does 
the original NFW profile. All of the halos presented by 
Navarro et al. (2004) are well described by the new pro- 
file. A key difference between this profile and the classic 
NFW or Moore et al. (1999) profiles is that there is no 
well-defined central cusp. The density reaches a finite 
central value rather than diverging. On observationally 
relevant scales (100 to 1000 pc), though, the density pro- 
files remain rather steep. Observations like ours with a 
resolution of tens or hundreds of parsecs would be ex- 
pected to find density cusps only marginally shallower 
than an NFW profile. Only on sub-parsec scales does 
the logarithmic slope of this profile reach even a = 0.5 
(Navarro et al. 2004). 

When the fits to Equation 3 for various halos are 
rescaled by the characteristic density {p-2) and radius 
(r_2) of each halo, they are all essentially identical, cor- 
roborating the hypothesis of a universal density profile. 
Diemand et al. (2004) confirm the ability of this func- 
tional form to fit their own independently simulated ha- 
los. However, Diemand et al. also show that a general- 
ized three-parameter NFW profile with a formal central 
cusp is able to fit the halos nearly as well. Note that our 
observations still probe well below the scales that are 
resolved by the simulations, and it is unknown which, if 
any, of the profiles motivated by the simulations provides 
the most accurate extrapolation to smaller radii. 

4.1.3. Comparison 

Despite our lack of knowledge about the precise func- 
tional form of simulated density profiles at small radii, 
we can still compare our observational results with the 
profiles that fit the simulations best. In a general sense, 
it is obvious merely from inspection of Figures 4 to 7 
(and Figure 10 in Simon et al. 2003) that the five ro- 
tation curves have rather different shapes. This visual 
impression is confirmed by the fits we performed in §3. 
If we use power laws to describe the density profiles, the 
mean slope is Odm = 0.73, with a dispersion of 0.44. By 
comparison, Diemand et al. (2004) find for the general- 
ized NFW profile fits that the central power law index is 
c^DM = 1.16±0.14. A more dramatic (although less intu- 
itive) illustration of the difference between the observed 
and simulated density profiles is provided by our fits with 
the new Navarro et al. (2004) density profile. Although 
this formula fits our data well, as shown in §3.3.2 and 
Table 6, we derive values of 77 ranging from 0.12 to 8.99. 
Even if we ignore the clearly absurd results for the flat 
density profile of NGC 2976, we find a mean value of 
rj — 0.73, with a dispersion of 0.78. From the simula- 
tions, Navarro et al. (2004) measure 77 = 0.172 ± 0.032 

Note that the highest-resolution simulations currently have a 
resolution limit of 500 — 1000 pc, so another factor of 5 improve- 
ment in the resolution is needed before the simulations reach the 
scales probed by observations. 
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and Diemand et al. (2004) arc in agreement, finding 
r) = 0.186 ± 0.037. Witli the exception of NGC 5963, 
our fits are seriously discrepant with these results. If we 
force 7] to lie within the range preferred for the simu- 
lated halos, the scale radius and scale density run away 
to unreasonable values (as can be seen in Table 6 for 
NGC 6689). 

Thus, there are two notable difFercnccs between our 
results and the most recent CDM simulations: 1) the 
halo-to-halo scatter is > 3 times larger than the simula- 
tions, and 2) the observed central slopes are on average 
shallower than the simulations. Both of these differences 
present challenges for future simulations. 

Although the shapes of the central density profiles we 
derive disagree with the theoretical predictions, the ac- 
tual values of the central densities are relatively con- 
sistent with the simiilation results. In Table 7 we give 
measured values for the dimensionless quantity Ay/2 in- 
troduced by Alam, Bullock, & Weinberg (2002) to pa- 
rameterize halo central densities. In a ACDM cosmology, 
the galaxies in our sample should have Ay/2 values be- 
tween 10^ and 3 x 10^ (assuming NFW density profiles), 
similar to the values we observe. The galaxies do have 
a tendency to lie somewhat below the predicted level, 
but the difference is within the la scatter expected for 
ACDM. 

4.2. The Significance of NGC 5963 

Despite these differences from the simulations, the den- 
sity profile of NGC 5963 appears to be in good agreement 
with the theoretical models. The rotation curve of this 
galaxy rises so rapidly than an NFW profile fits the data 
very well. A power law with a slope steeper than NFW 
provides an even better fit. To our knowledge, NGC 5963 
is the only low-mass disk galaxy for which a CDM-like 
central density cusp is confirmed and shallow density pro- 
files are ruled out. 

Since NGC 5963 represents the exception rather than 
the rule, the question is, are all galaxies formed with 
cuspy density profiles and most lose them, or did an un- 
usual event during the evolution of NGC 5963 caused its 
density profile to become so steep? The key to answering 
this question is identiiying what about NGC 5963 makes 
it unique. 

A number of galaxy properties from our sample 
are summarized in Table 7. Among these galaxies, 
NGC 5963 has the highest central surface brightness, de- 
spite its very low surface brightness outer disk. It also has 
the most unusual surface brightness profile (see §3.5.3), 
although the significance of this is not clear. NGC 5963 
is at the brighter end of the luminosity range of our tar- 
gets, but the total span between the five galaxies is only 
~ 1 magnitude. It does have the largest rotation veloc- 
ity, and therefore the highest mass, among our sample. 
Other than mass, the most obvious distinction between 
NGC 5963 and the other galaxies is the sharp transi- 
tion between the inner and outer surface brightness pro- 
files. This causes the derived stellar disk rotation curve 
to peak at small radii (r = 20") and then drop steeply, 
making the outer parts of the galaxy highly dark matter 
dominated (see Figure 6b). Although this photometric 
structure may be related to the steep density profile of 
NGC 5963, we do not know its physical origin. If the 



galaxy really does contain a pseudobulge, then the pro- 
cesses that have funneled baryons toward the center of 
the galaxy could also have caused the dark matter to 
become more centrally concentrated. 

It is also noteworthy that high resolution two- 
dimensional velocity fields and deep multicolor imaging 
are not necessary to recognize the unique dark matter 
density profile of NGC 5963. Using lower resolution H I 
data, a long-slit spectrum of modest resolution, and pho- 
tographic U and R band imaging more than 15 years ago, 
Bosma et al. (1988) already concluded that this galaxy 
probably contains "an unusually centrally concentrated 
and massive halo" . Since high-resolution CDM simula- 
tions of galaxy structure had not yet been done Bosma 
et al. could not generalize this result into any broader 
implications for dark matter or galaxy formation. 

4.3. NFW and Pseudo-Isothermal Fits to Rotation 

Curves 

Our study is not alone in finding that pseudo- 
isothermal profiles fit many rotation curves better than 
do NFW profiles. A number of previous authors have 
reached similar conclusions (e.g., de Blok et al. 2001a; 
Borriello & Salucci 2001; dc Blok, McGaugh, & Rubin 
2001b; de Blok & Bosma 2002; Salucci, Walter, & Bor- 
riello 2003; Simon et al. 2003). We emphasize, however, 
that the observation that pseudo-isothermal fits produce 
lower reduced values than NFW or Moore profiles 
does not demonstrate the presence of a constant-density 
core. We showed in §3.3 that in many cases power law 
density profiles with intermediate (0 < a < 1) slopes 
provide comparable fits. A modified pseudo- isothermal 
profile with an a = 0.5 central cusp (Equation 2) also fits 
these rotation curves as well as a true pseudo-isothermal 
profile with a const ant- density core does. Since current 
theories of galaxy formation do not lead naturally to ei- 
ther cores or pseudo-isothermal density profiles, a prefer- 
ence for constant-density cores over shallow density cusps 
is not justified. Whether these galaxies actually contain 
a central region of constant density or if they have shal- 
low cusps cannot be determined from the present data 
despite the high resolution of our measurements. 

We do agree with past studies that for most galaxies 
the NFW form does not fit the inner density profile very 
effectively (see Figures 4b- 7b). When rg and r2oo are 
constrained to stay roughly in the range expected from 
simulations, the values for NFW fits are generally 
rather high. We also point out, however, that the NFW 
fit parameters are quite poorly constrained unless the 
observations extend well into the flat part of the rotation 
curve. If there are only a few data points on the flat part 
of the rotation curve (especially when these data points 
have the largest uncertainties, as is often the case), they 
do not provide significant leverage on the fit. In such 
cases, the NFW scale radius and virial radius (or any 
other pair of parameters that can be used to describe an 
NFW profile) become completely degenerate. In order 
to obtain accurate estimates of the NFW fit parameters, 
our results suggest that > 50 % of the observed part of 
a galaxy must have a nearly flat rotation curve. 

4.4. Implications for CDM 

Combining this study with our previous analysis of 
NGC 2976, we have shown that the five galaxies in our 
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TABLE 7 

Are Galaxy Parameters Correlated With Density Profile Slope? 



Galaxy 


Distance 
[Mpc] 


Mi 
[mag] 


[mag aJTcsec"^] 


V -I 
[mag] 




M,/Lr 


i 

[deg] 


T/ C 

^ max , halo 

[km s-l] 


[pc] 






NGC 2976^ 


3.45 


-18.6 


19.73 


0.87 


3.7 X 10^ 


0.53 


61.4 


74 


900 


7.0 X 10"' 


0.01 


NGC 4605 


4.26 


-18.9 


19.36 


0.75 


4.7 X 10-' 


1.01 


71.5 


74 


560 


1.7 X 10^' 


0.78 


NGC 5949 


14.0 


-19.8 


19.75 


0.95 


1.1 X 10^0 


1.64 


64.6 


74 


880 


6.8 X 10"' 


0.88 


NGC 5963 


13.0 


-19.1 


18.10 


0.85 


1.4 X 10^° 


1.24 


48.4 


114 


660 


2.9 X 10'5 


1.20 


NGC 6689 


11.0 


-19.0 


21.25 




1.5 X IQio 


1.968 


76.0 


94 


1330 


4.8 X 106 


0.79 



''Extrapolated central surface brightness of the disk. 

''We calculate the dynamical mass as M^y„ = V^^^Rmax/G, where Vmax is the maximum observed rotation velocity and Rmax is the 
largest radius our observations reach. Since the galaxies extend to much larger radii, these numbers clearly represent lower limits to the 

actual masses of each galaxy. 
'^Thc maximum rotation velocity of the dark matter halo after removing the stellar disk. 

"^The radius at which the dark matter rotation curve reaches half of Vj^ax haloi ^ defined by Alam et al. (2002). 
"The mean dark matter density within Ry/2, in units of the critical density. 
^Using data from Simon et al. (2003). 

SThis mass-to-light ratio is calculated for the Sloan r' band, not the Kron-Cousins R band. 



sample have density profiles with a wide range of central 
slopes. If we fit the data with power laws, we find that 
three galaxies have central cusps of cxdm ~ 0.8, while 
the other two galaxies have very different central slopes 
(aoM = 0.01 and a^M = 1-20). Of the outliers, NGC 2976 
is a satellite of M 81, and thus may have been subject 
to tidal stripping. There are suggestions in the literature 
that tidal stripping acts to remove density cusps (Stoehr 
ct al. 2002; Hayashi ct al. 2003), but the most recent 
simulations find that the central density slope is not al- 
tered by tidal effects (Kazantzidis et al. 2004b). It is not 
clear why NGC 5963 differs so strongly from the other 
galaxies. These fits demonstrate that while galaxies with 
steep central density cusps do exist, they arc not shared 
by all galaxies. Furthermore, no single value of the cen- 
tral slope can describe all five of the rotation curves in 
our sample. 

On the other hand, if we consider the pseudo- 
isothermal fit results (as most previous observational 
studies have done), we find that good fits can be ob- 
tained for four of the five galaxies (sec Table 4). It is 
interesting to note that these galaxies all have central 
densities within a factor of seven of each other. Only 
for NGC 5963 is the quality of the pseudo-isothermal fit 
sufficiently poor that a constant-density core can be ex- 
cluded. The rotation curves of the other four galaxies are 
consistent with cores, but all except NGC 2976 are also 
consistent with cusps. The implication of these findings 
is that some previous studies may have overestimated the 
disagreement between their data and the CDM simula- 
tions. 

Based on these results, we reach the following conclu- 
sions: 

• First, the observed variation in density profile 

slopes from halo to halo is much larger than ex- 
pected from the simulations. We also find that 
none of the commonly used density profile func- 
tional forms can describe all five galaxies. 

• Second, most galaxies have density profiles signif- 
icantly shallower than the q;dm = 1 central cusps 
preferred by the simulations. Nevertheless, the dif- 



ference between the central slopes we measure and 
the theoretical values is much smaller than sug- 
gested by most previous studies. 

• And third, some galaxies do contain cvdm 

cusps with an overall density profile that follows 
the NFW form, hinting that it may still be possi- 
ble to reconcile the results of the pure dark matter 
cosmological simulations with observations. 

Do these differences between the observations and the 
simulations indicate a fundamental problem with the 
CDM paradigm? Probably not. A great many plausible 
theoretical explanations have been proposed in attempts 
to understand this problem, and a few of the recent ideas 
are mentioned below. Taylor & Navarro (2001) suggested 
based on phase-space arguments that CDM density pro- 
files should have a = 0.75 cusps, very similar to the 
average we measure. Ricotti (2003) found that the dark 
matter halos of simulated dwarf galaxies have shallower 
cusps than those of massive galaxies. In a similar ex- 
periment, however, Colin et al. (2004) reached the oppo- 
site conclusion, so the degree to which dwarf galaxy den- 
sity profiles are expected to match those of large galaxies 
has not yet been settled. Boylan-Kolchin & Ma (2004) 
showed that major mergers leave cuspy remnants unless 
both precursors had cores, so at least in the case of pure 
dark matter halos, mergers cannot destroy cusps. Bary- 
onic processes, however, probably can (e.g., Weinberg & 
Katz 2002). Even the dark matter itself may be able 
to flatten cuspy density profiles; Ma & Boylan-Kolchin 
(2004) argued that energy deposition by merging dark 
matter substructures can puff up central cusps into shal- 
lower density profiles. This effect might have escaped 
notice in previous simulations due to insufficient resolu- 
tion or overmerging (unrealistic destruction of subhalos 
in dense regions). Finally, Hayashi et al. (2004) have 
shown that if dark matter halos are significantly triaxial 
(see §5) , for some viewing angles the derived density pro- 
files can appear much shallower than the actual density 
profiles. 

Since we have shown that cuspy density profiles are 
consistent with a large fraction of our sample, and in light 
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of the variety of potential effects not currently accounted 
for in the simulations that could change the predicted 
density profiles, there is little reason to conclude at this 
point that the density profile controversy represents a 
crisis for CDM. 

5. THE EFFECTS OF HALO TRIAXIALITY 

Almost every previous observational study of density 
profiles and rotation curves has assumed that dark mat- 
ter halos arc spherical. CDM simulations, however, sug- 
gest that dark matter halos are triaxial (Dubinski & Carl- 
berg 1991; Warren ct al. 1992; Cole & Lacey 1996; Jing & 
Suto 2002). Unfortunately, very few observations of indi- 
vidual galaxies are available to test this prediction. Sev- 
eral polar-ring galaxies have been studied, with results 
ranging from axis ratios as small as 0.3-0.4 up to axis ra- 
tios of 1 (spherical halos) (Schweizer, Whitmore, & Ru- 
bin 1983; Whitmore, McElroy, & Schweizer 1987; Sackett 
& Sparke 1990; Sackett et al. 1994; Combes & Arnaboldi 
1996; lodice ct al. 2003). A variety of techniques sug- 
gest that the gravitational potential of the dark matter 
halo of the Milky Way is close to spherical (Kuijken & 
Tremaine 1994; Oiling & Merrifield 2000). Despite the 
scarcity of reliable constraints on individual halo elliptici- 
ties, statistical arguments provide a way to determine the 
average shape of galaxy halos. Franx & de Zeeuw (1992) 
show that the observed scatter in the Tully-Fishcr (TF) 
relation places a strong upper limit on the allowed el- 
lipticities of spiral galaxy disks. Even if all of the TF 
scatter is caused by elongated disks, the mean disk el- 
lipticity is required to be less than 0.1. More likely, the 
scatter is a result of a combination of effects, indicating 
that on average disk ellipticities are 0.05 or smaller. In 
this section we consider the effects triaxial halos could 
have on the velocity fields we observe. 

5.1. Halo Oblateness 

Galaxy disks are expected to be oriented such that the 
disk lies in the plane described by the major and inter- 
mediate axes of the halo, with the z-direction pointing 
along the halo minor axis (Sharma & Steinmetz 2004; 
Bailin & Steinmetz 2004). We first consider the effect 
of this flattening of the halo, under the assumption that 
halo shape in the disk plane is circular. 

The rotation velocity of a disk in such an oblate 
spheroidal halo is given by Equation 2-91 in Binney & 
Tremaine (1987): 

2 /n\ A ^ 7 /"^ p{ni)nrdm 

where p{m) is the density profile of the halo in the 
spheroidal coordinate m (defined by rn? = r'^ + /{I — 
e^)) and e is the eccentricity of the spheroid. For any 
density profile, this integral can be evaluated numeri- 
cally to give the rotation curve for various degrees of 
oblateness. We performed this calculation for spheroidal 
power laws {p{m) = /9o(™/'Tio)~") and NFW profiles 
{p{m) = pcScim./m.s)^^ (1 + (m/m.s)^)^^). The change 
from a spherical halo to an oblate one could affect both 
the shape of the rotation curve and its amplitude. We 
find that as the halo becomes fiatter (e — > 1) the ampli- 
tude of the rotation curve increases (as expected, since 
the mass becomes concentrated closer to the disk). We 



also find that the shape of the rotation curve is com- 
pletely independent of the halo fiattening for a power 
law; the same power law index for the rotation curve 
(and density profile) is derived for any value of e. For 
an NFW density profile the shape of the rotation curve 
changes very subtly with eccentricity (the peak of the 
rotation curve shifts to smaller radii as the halo becomes 
flatter), but the inner slope of the rotation curve is es- 
sentially unaffected. 

5.2. Disk Ellipticity 

In the previous subsection, we showed that the flatten- 
ing of the halo does not alter the observed density profile 
for a circular disk. We now remove the assumption that 
the disk is circular and study the effects of disk ellipticity 

on the observed velocity field. 

5.2.1. Are Noncircular Motions Common in Disk 
Galaxies? 

Four of the five galaxies in our sample show evidence 
for a radial component to their velocity fields. In two of 
these, the radial motions are detected very strongly; for 
the other two the significance of the radial term is lower, 
although it is still confirmed at 95 % confidence. Other 
recent studies have also begun to find significant numbers 
of galaxies with noncircular motions (Schoenmakers et 
al. 1997; Swaters et al. 2003; Coccato et al. 2004; Blais- 
Ouellette ct al. 2004; Wong, Bhtz, & Bosma 2004). A 
large majority of the galaxies studied by these authors 
(and us) are not barred. However, of the galaxies in 
which observations could have detected radial motions, 
nearly all indeed appear to contain them. Are radial 
components to the velocity fields ubiquitous in late-type 
spiral galaxies? What is the cause of these motions? 

5.2.2. Measuring Disk Ellipticity With Noncircular 

Motions 

Since there is no strong evidence for bars in the four 
galaxies where we detect noncircular motions, other 
sources of these motions should be considered. First, 
we note that the observed noncircular motions are domi- 
nated by the radial component. The simplest possibility 
is that the galaxies contain strong radial fiows directed 
towards their centers. Given the observed magnitude of 
the radial motions at a radius of ~ 1 kpc and an assumed 
volume density for the gas (1 H atom cm~^), if the radial 
motions actually represent a net inflow of gas we conser- 
vatively estimate that all of the gas would accumulate 
within a ~ 1 kpc radius of the galaxy centers in 1-3 Gyr. 
The star formation rate over the same region is at least 
an order of magnitude too small to consume th(^ inflow- 
ing gas. We therefore conclude that this interpretation is 
not viable. The most intriguing remaining explanation 
is that the gas is moving on elliptical orbits, which could 
result from the influence of a triaxial dark matter halo. 

In a disk galaxy with a triaxial dark matter halo, the 
potential in the plane of the disk will in general be el- 
liptical. The closed orbits in such a potential are also 
ellipses, which means that the observed line-of-sight ve- 
locities of an elliptical disk will differ from pure circular 
rotation. This problem has been considered in detail by 
Binney (1978), Teuben (1991), Franx, van Gorkom, & de 
Zeeuw (1994), Schoenmakers et al. (1997), and Schoen- 
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makers (1998). Franx ct al. (1994) showed that an cUip- 
tical potential induces components in the velocity field 
proportional to the ellipticity and the angle between the 
long axis of the ellipse and the observer's line of sight. 
Adopting their formalism (Pranx et al. 1994, Equation 
AlO), we have the following relations: 



Cl 



1 - (- - a)eR cos 2(j)obs 



4(1 - g2 



a tR sin 2(/)o6s 



fcSini (6) 
VcSini, (7) 



where ci and si are the cos9 (rotation) and svnO (radial) 
components of the tilted-ring model, respectively. The 
hat symbols indicate that these quantities are derived un- 
der the assumption of circular orbits, and will therefore 
differ systematically from the true values if the potential 
is elliptical. These formulae are valid when the rotation 
curve, v^(r), can be represented as a power-law with in- 
dex (3, which we have shown is a reasonable assumption 
for these galaxies. (3 is related to a by a = ^/3/(l + /3). 
q is the axis ratio of the galaxy, which we take from the 
photometry, (pobs is the angle between the long axis of 
the potential and the line of sight, and en is the elliptic- 
ity of the orbit. The relationship between the ellipticity 
of the potential, Cpoti and that of the orbit, e^, is given 
by Pranx et al. (1994): 



^pot 



1-/3 
1 + /3 



(8) 



In order for us to make the simplifying assumption that 
the apparent rotation velocities are equal to the true ro- 
tation velocities (i.e., that the orbits are nearly circular), 
the following condition must be met: 



(I - a)eflcos2^o6s «; 1. 



(9) 



In this case. Equation 6 reduces to ci = Vcnini. Sub- 
stituting this expression into Equation 7 and solving for 
the potential ellipticity yields 



e_R sin 2(pobs 



1 
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Cl 



(10) 



Thus, if the observed noncircular motions are indeed 
due to an elliptical potential in the plane of the disk, 
the amplitudes of the circular and radial motions put 
a direct constraint on the ellipticity. Note that since 
sin 2(j)obs < 1 this constraint is actually a lower limit on 
e_R. Using the known values of q and a, and the re- 
sults of our tilted-ring models for ci and si, we calcu- 
late Eij sin 2(^obj, as a function of radius for each galaxy. 
The results are shown in Figure 9. NGC 4605, with its 
very strongly detected radial motions, has a large and 
nearly constant ellipticity that deviates significantly from 
zero beginning at a radius of 20". The mean value of 
en sin 2^obs over the entire galaxy is 0.175 ± 0.016, mak- 
ing the disk of this galaxy substantially elliptical. From 
Equation 8, we see that epot = O.lSen, so the elliptic- 
ity of the potential is at least 0.03. eiisiii2(j)obs is de- 
tected at the ~ So- level in NGC 5949 and NGC 5963. 
We find that €Rsm2(l)obs = 0.043 ± 0.014 for NGC 5949 
and eRsm24>obs = 0.060 ± 0.020 for NGC 5963. We do 
not detect any evidence for elliptical orbits in NGC 6689: 
CR sin 2(f>obs = 0.007±0.014. NGC 2976 is the only galaxy 



in which the ellipticity is not constant with radius. This 
could suggest that either the halo structure is more com- 
plicated in this galaxy (e.g., the axis ratios of the halo 
change with radius), or that its radial motions may be 
caused by something other than triaxiality. If the halo 
of NGC 2976 is triaxial, the mean ellipticity of the disk 
orbits is eRsm2(t>obs = 0.116 ±0.013. 

Given the determination of eR sin 2(j)obs y we can now go 
back and confirm that the assumption stated in Equation 
9 (that the difference between the apparent and true ro- 
tation velocities is small) was justified. As long as the 
ellipticity is modest {er < 0.2, which requires that the 
viewing angle be larger than ^ 9° for most of these galax- 
ies), this condition will be satisfied. Note, however, that 
for NGC 4605 cr may be larger than this value. 

These calculations show that, if the observed radial 
motions are indeed due to halo triaxiality, the lower limit 
on the mean potential ellipticity is of order a few per- 
cent. Since the scatter in the TF relation indicates that 
the average ellipticity of the potential has a strict upper 
limit of 0.1 (Pranx & de Zeeuw 1992), the ellipticities 
are now constrained within a relatively narrow range. 
CDM simulations predict that the dark matter halos of 
galaxies similar in mass to our targets should have axis 
ratios of c/a = 0.47 and h/a = 0.62 (Jing & Suto 2002), 
which would suggest larger ellipticities than we measure. 
Bailin & Steinmetz (2004) find somewhat less triaxial ha- 
los with c/a = 0.65 and b/a = 0.83 for halos with masses 
of ^ 3 X IQ^^Mq. Adding gas cooling to the simulations, 
however, may make the inner regions of the halo much 
more spherical (Kazantzidis et al. 2004). With obser- 
vations of a larger sample of galaxies, it may be possi- 
ble to obtain robust constraints on the three-dimensional 
shapes of galaxy halos, which would then provide another 
strong test of CDM simulations on small scales. 

6. SUMMARY AND CONCLUSIONS 

We have used two-dimensional CO and Ha velocity 
fields, sampled at high spatial resolution and high spec- 
tral resolution, to constrain the dark matter density pro- 
files of four nearby, low-mass, late-type galaxies. We 
obtained rotation curves from the data cubes by con- 
structing tilted-ring models, and found that three of the 
four galaxies contain radial motions in addition to ro- 
tation. Combining these data with the observations of 
NGC 2976 presented by Simon et al. (2003) gave us a 
sample of five galaxies to study. 

For each galaxy we constructed a model of the stel- 
lar disk and investigated the density profile of the dark 
matter halo under varying assumptions about the stel- 
lar mass-to-light ratio. In most cases, over the range 
of plausible mass-to-light ratios, the dark matter den- 
sity profiles do not change significantly. For the mass- 
to- light ratios indicated by the galaxy colors, we found 
that the five galaxies span a range of central density pro- 
file slopes from = to Odm = 1-20. The mean 
slope is a-om = 0.73, with a dispersion of 0.44. Neither 
of the standard density profile functional forms (pseudo- 
isothermal and NEW) is able to adequately fit all five 
galaxies. The scatter in slope that we observe is three 
times larger than that seen in the simulations, and the 
mean slope is smaller than predicted. We do note, how- 
ever, that NGC 5963 is the first low-mass disk galaxy 
in which a cuspy density profile of the predicted form is 
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Fig. 9.— €i{sin2(^<,6s as a function of radius for (from top to bottom) NGC 2976, NGC 4605, NGC 5949, NGC 5963, and NGC 6689. 
The shaded gray bands across each panel represent the weighted mean value of e^; sin2^o(,3. The points at large radii without error bars 
correspond to rings in which the radial velocity was held fixed because the ring contained too few data points or the angular coverage of 
the data was not sufficient to constrain the radial motions. These points do not contribute to the weighted average and should not be 
considered to contain significant information about the orbital ellipticity. Only in NGC 6689 does the gas appear to be traveling on circular 
orbits, although it is also possible that cp^bs ^ 0° oi^ 90° for this galaxy. The ellipticity of the potential is a factor of a few smaller than 
the ellipticity of the orbits that is plotted here (see Equation 8). Note that these values are lower limits on the true ellipticity due to the 
unknown viewing angle. 



confirmed, while shallower alternative profile shapes are 
ruled out. 

Although these results still indicate significant dis- 
agreements between CDM simulations and the properties 
of real galaxies, there are also some differences between 
the findings of this paper and those of previous obser- 
vational analyses. Many authors have concluded that 
density profiles with constant-density cores fit the rota- 
tion curves of dwarf and LSB galaxies better than cuspy 
density profiles do (e.g., de Blok ct al. 2001a; Borriello 
& Salucci 2001; de Blok et al. 2001b; de Blok & Bosma 



2002; Salucci et al. 2003; Simon et al. 2003). Typi- 
cal results from these studies are that the distribution 

of density profile slopes is a = 0.2±0.2 (de Blok et 
al. 2001a; de Blok & Bosma 2002; de Blok et al. 2003). 
Some other studies agreed that cored profiles generally 
fit the data better, but argued that this effect is pri- 
marily the result of systematic uncertainties in the ob- 
servations (van den Bosch et al. 2000; van den Bosch 
& Swaters 2001; Swatcrs et al. 2003). When we fit 
power laws to the rotation curves we find significantly 
steeper density profile slopes on average than previous 
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authors did, although both pseudo-isothermal density 
profiles (constant-density cores) and relatively cuspy pro- 
files are consistent with our data. Despite our disagree- 
ment on the value of the mean slope, we agree with other 
recent studies that found that dark matter density pro- 
files span a wide range of slopes (e.g., de Blok et al. 2003; 
Swaters et al. 2003). 

What is responsible for the steeper slope that we find? 
One important element of our study is that the two- 
dimensional velocity fields, high spatial resolution, and 
target selection we employed significantly reduce the im- 
pact of systematic uncertainties on our results (Simon 
et al. 2003). As an illustration of this point, van den 
Bosch & Swaters (2001) showed that to obtain tight con- 
straints on density profile slopes, at least one of the three 
following conditions must be met: 1) the uncertainties on 
the rotation curve must be < 0.2 km s^^, 2) the rotation 
curve must extend out beyond 20% of the virial radius, 
and 3) the rotation curve must extend inward to radii less 
than 2% of the virial radius. While we do not reach 
the velocity precision of criterion (1), we approach it as 
closely as feasible, with typical systematic uncertainties 
as small as ~ 2 km s^^. Although higher velocity resolu- 
tion spectroscopy is possible, our velocity measurements 
have already run into the limit imposed by the small-scale 
random motions of ~ 5 km s~^ that most galaxies appear 
to contain. Without H I observations at large radii, we 
also fail to meet criterion (2), but our high spatial resolu- 
tion does allow us to satisfy criterion (3). Galaxies with 
masses similar to those studied here may have virial radii 
of 50 - 75 kpc (Bullock et al. 2001; Navarro et al. 2004), so 
even for our most distant targets our ^ 250 pc resolution 
probes well within 2% of the virial radius. In addition, we 
explored functional forms for the density profiles other 
than the standard NFW/pseudo-isothermal dichotomy, 
thereby demonstrating that good pseudo-isothermal fits 
in general do not rule out cuspy density profiles. Given 
these new results, we conclude that while the discrepan- 
cies between CDM and galaxy rotation curves may not be 
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as severe as previously thought, improved simulations are 
needed in order to determine whether the theoretical and 
observational findings can be brought into agreement. 

Finally, we considered the origin of the noncircular mo- 
tions that are now being detected regularly in disk galax- 
ies. A number of previous authors have showed that a 
triaxial dark matter halo can produce exactly the ob- 
served effect on the gas orbits. Under the assumption 
that halo triaxiality is causing the noncircular motions, 
we measure lower limits to the orbital ellipticities of 0.02- 
0.17. These correspond to lower limits on the ellipticities 
in the potential between 0.005 and 0.03. 
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